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ABSTRACT 

An  expository  and  illustrated  treatment  of  the  basic  B-spline  theory  as  derived  from 
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SIGNIFICANCE  AND  EXPLANATION 


This  report  contains  the  lecture  notes  for  the  (irsl  of  four  lectures  which  comprise 
the  course  entitlecT“Tho  extension  of  B-spline  curve  algorithms  to  surfaces*'  given  at  SIG- 
GRAPH’86.  It  is  an  elaboration  am!  extension  of  the  MRC  report  ^2896  by  de  Boor 
and  Hollig.  in  which  the  basic  B-spline  theory  is  developed  from  the  recurrence  relation 
rather  than  the  original  definition  in  terms  of  divided  differences  of  the  truncated  power. 
This  avoids  what,  to  the  people  in  CAGD.  amounts  to  a  detour  through  the  theory  of 
divided  differences.  Somewhat  surprisingly,  the  resulting  development  is  no  longer  than 
the  standard  one.  and  in  some  respects  seems  even  more  direct.  It  does  bring  to  the  fore 
the  dual  functionals  and  stresses  the  point  that  B-splines  are  best  treated  in  terms  of  their 
linear  span. 
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The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive  summary  lies 
with  MRC.’.  and  not  with  the  author  of  this  report. 
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0.  Introduction 

These  lecture  notes  review  those  basic  lacts  about  (univariate)  B-splines  which  are 
of  interest  in  CAGD.  The  intent  is  to  give  a  self-contained  and  complete  development  of 
the  material  in  as  simple  and  direct  a  way  as  possible.  For  this  reason,  the  B-splines  are 
defined  via  the  recurrence  relations,  thus  avoiding  the  discussion  of  divided  differences 
which  the  traditional  definition  of  a  B-spline  as  a  divided  difference  of  a  truncated  power 
function  requires.  As  this  lecture  is  intended  to  show,  this  does  not  force  more  elaborate 
derivations  than  are  available  to  those  who  feel  at  ease  with  divided  differences.  It  does 
force  a  change  in  the  order  in  which  facts  are  derived  and  brings  more  prominence  to  such 
things  as  Marsden’s  Identity  or  the  Dual  Functionals  than  they  currently  have  in  CAGD. 

In  addition,  it  highlights  the  following  point:  The  consideration  of  a  single  B-spline 
is  not  very  fruitful  when  proving  facts  about  B-splines,  even  if  these  facts  (such  as  the 
smoothness  of  a  B-spline)  can  be  stated  in  terms  of  just  one  B-spline.  Rather,  simple 
arguments  and  real  underst  anding  of  B-splines  are  available  only  if  one  is  willing  to  consider 
all  the  B-splines  of  a  given  order  for  a  given  knot  sequence.  Thus  it  focuses  attention  on 
splines,  i.e.,  on  the  linear  combination  of  B-splines. 

The  lecture  deals  with  splines  for  an  arbitrary  knot  sequence  and  does  rarely  become 
more  specific.  In  particular,  the  B(ernstein-Bezier)-net  for  a  piecewise  polynomial,  though 
a  (very)  special  case  of  a  representation  by  B-splines.  gets  much  less  attention  than  it 
deserves,  given  its  immense  useiulness  in  CAGD  (and  spline  theory).  But  the  third  lecture 
takes  up  this  topic. 

The  lecture  deals  only  with  spline  functions.  There  is  an  immediate  extension  to 
spline  curves:  Allow  the  coefficients,  be  they  B-spline  coefficients  or  coefficients  in  some 
polynomial  form,  to  be  points  in  IR*  or  IR3.  But  this  misses  the  much  richer  structure 
for  spline  curves  available  because  even  discontinuous  parametrizations  may  describe  a 
smooth  curve.  This  topic  of  geometric  continuity  is  discussed  in  deta.l  in  the  fourth 
lecture. 

The  lecture  notes  are  solidly  based  on  1BH86-  which  covers  more  or  less  the  same 
material,  in  a  Ies^  elaborate  way  and  without  any  figures,  in  just  seven  pages. 

The  relevant  literature  on  (univariate)  B-splines  up  to  about  1975  is  summarized  in 
B76  which  also  contains  hints  of  the  most  exciting  developments  concerning  B-splines 
since  then:  knot  insertion  and  the  multivariate  B-splines.  These  are  covered  in  the  second 
lecture,  but  knot  insertion  is  already  put  to  good  use  in  the  last  part  of  this  lecture.  The 
two  books  on  splines.  B78  and  SchuSl  .  which  have  appeared  since  1975.  cover  B-splines 
in  the  traditional  way.  As  presentations  of  splines  from  the  CAGD  point  of  view,  the 
survey  article  BFK84  and  the  “Killer  B  s"  BBB85.869  are  particularly  recommended. 
Sponsored  bv  the  United  States  Army  under  Contract  No.  DA  AC.  29- 80-004 1 . 
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1.  B-splines  defined 

We  start  with  a  partition  or  knot  sequence,  i.e.,  a  nondecreasing  sequence  t  := 
(t,).  The  B-splines  of  order  1  for  this  knot  sequence  are  the  characteristic  functions  of 
this  partition,  i.e.,  the  functions 


Bii(t)  :=X.(t)  :=  {*’  ,f[> 

'  '  '  '  \  0,  otherwise. 


(1-1) 


Note  that  all  these  functions  have  been  chosen  here  to  be  right -continuous.  Other  choices 
could  have  been  made  with  equal  justification.  The  only  constraint  is  that  these  B-splines 
should  form  a  partition  of  unity,  i.e., 


Eft-  (f)  =  1,  for  all  t. 


(1.2) 


In  particular, 


<,  =  *,+,  =*  Bti  =  X,  =  0.  (1.3) 

From  these  first-order  B-splines,  we  obtain  higher-order  B-splines  by  recurrrence: 


with 


Vi+i,k)Bi+i,k-i 

(1.4a) 

if  U  ^  1 

otherwise. 

(1.46) 

Thus,  the  second-order  B-spiine  is  given  by 

B,2  =  u>l2X,  +  (1  —  Wt+i.jJX,*!, 


(1.5) 


and  so  consists,  in  general,  of  two  nontrivial  linear  pieces  which  join  continuously  to  form 
a  piecewise  linear  function  which  vanishes  outside  the  interval  [f,-,  f:+2(.  For  this  reason, 
some  call  Bx2  a  linear  B-spline.  If.  e.g.,  t,  =  +  j  (hence  X,  =  0).  but  still  <  f,.2. 

then  Bx 2  consists  of  just  one  nontrivial  piece  and  fails  to  be  continuous  at  the  double 
knot  t,  —  j,  as  is  shown  in  Fig.  1.1. 
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The  third-order  B-spline  is  given  by 

Bt3  =  *  (l  -  3)^1 +  1.2 

—  ^'t3+'i2^t  "*■  (Wts(l  -  ^’t  +  1,2)  +  (1  “  (1-6) 

(i  "  1 ,3 )  f  i  ~  W’,^2.2)Xt-,2 

This  shows  that,  in  general,  B,3  consists  of  3  (nontrivial)  quadratic  pieces,  and,  to  judge 
from  the  Fig.  1.2,  these  seem  to  join  smoothly  at  the  knots  to  form  a  C 1  piecewise  quadratic 
function  which  vanishes  outside  the  interval  lU,t\  +  3[-  Coincidences  among  the  knots 
,t, >3  would  change  this.  If,  e.g.,  t,  —  t,-+ 1  —  <t+2  (hence  X,  --  X,+  i  =  0),  then 
B,a  consists  of  just  one  nontrivial  piece,  fails  to  be  even  continuous  at  the  triple  knot  t,. 
but  is  still  C1  at  the  simple  knot  lt  r 3,  as  is  shown  in  Fig.  1.2. 


After  k  -  1  steps  of  the  recurrence,  we  obtain  B,^  in  the  form 

t  -*-k  l 

Bik  -  ^*X,.  (l-”i 

)  -■ 

with  each  b,^  a  polynomial  of  degree  <  k  since  it  is  the  sum  of  products  of  k  1  linear 
polynomials. 

From  this,  we  infer  that  B,*  is  a  piecewise  polynomial  of  degree  k  which  vanishes 
outside  the  interval  and  has  possible  breakpoints  t,,  •••,/,.  *.  In  particular.  B,* 
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is  just  the  zero  function  in  case  t,  -  Also,  by  induction,  Btk  is  positive  on  the  open 
interval  \tt,tl^k.,  since  both  u:,k  and  ]  —  are  positive  there. 


Figure  1.4  The  two  weight  functions  in  (1.4a)  are  positive  on  }tt,tt  +  k'-- 
supp  Blk. 

Further,  we  see  that  Blk  is  completely  determined  by  the  k- f  1  knots  <, - For 

this  reason,  the  notation 

B(-\tu...,tt+k)  :=  Blk  (1.8) 

is  sometimes  used.  Other  notations  in  use  include 

A\*  :=  Btk  and  Mlk  :=  [k/(U +*  -  tt))Blk.  (1.9) 

The  many  other  properties  of  B-splines  are  derived  most  easily  by  considering  not 
just  one  B-spline  but  the  linear  span  of  all  B-splines  of  a  given  order  k  for  a  given  knot 
sequence  t.  This  brings  us  to  splines. 


2.  Splines  defined 

A  spline  of  order  k  with  knot  sequence  t  is,  by  definition,  a  linear  combination 
of  the  B-splines  Blk  associated  with  that  knot  sequence.  We  denote  by 


•>k 


,t  Blka,  :  a,  ~  IR} 


(2.1; 


the  collection  of  all  such  splines. 

We  have  left  open  so  far  the  precise  nature  of  the  knot  sequence  t.  other  than  to 
specifv  that  it  be  a  nondecreasing  real  sequence.  In  any  practical  situation,  t  is  necessarily 
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a  finite  sequence.  But,  since  on  any  nontrivial  interval  +  i  at  most  k  of  the  Btk 
are  nonzero,  viz.  B}  ~k+i,k,  ■  •  • ,  Bjk,  it  doesn’t  really  matter  whether  t  is  finite,  infinite, 
or  even  bi-infinite;  the  sum  in  (2.1)  always  makes  point  wise  sense,  since,  on  any  interval 
<tj,tjJt  i j,  at  most  k  summands  are  not  zero. 


■  j 


•  J 
r’*« 
A* 


Figure  2.1  The  k  B-splines  whose  support  contains  ,t;, t}n.x 
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We  will  pay  special  attention  to  the  following  two  “extreme”  knot  sequences,  the 
sequence 

TL  :=  (...,  -2,  -1,0, 1,2,...) 


and  the  sequence 


LB  :=  (. . .  ,0.0,0. 1,1.1....) 


A  spline  associated  with  the  knot  sequence  2Z  is  called  cardinal  splines.  This  term 
Weis  chosen  by  Schoenberg  Scho(i9  because  of  a  connection  to  Whittaker's  Cardinal  Series. 
This  is  not  to  be  confused  with  its  use  in  earlier  spline  literature  where  it  refers  to  a  spline 
which  vanishes  at  all  points  in  a  given  sequence  except  for  one  at  which  it  takes  the  value 
1.  The  latter  splines,  though  of  great  interest  in  spline  interpolation,  do  not  interest  us 
here. 

Because  of  the  uniformity  of  the  knot  sequence  t  -  TI..  formulae  involving  cardinal 
B-spiines  are  often  much  simpler  than  corresponding  formulae  for  general  B-splines.  To 
begin  with,  all  cardinal  B-sj>l  tries  (of  a  given  order  I  are  translates  of  one  another.  With 
the  natural  indexing  t ,  for  the  entries  of  the  uniform  knot  sequence  t  =  TL.  we 

have 


\k  -vA-(-  -  »). 

(2-2) 

B.k  B(  •() . k). 

(2.3) 
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The  recurrence  relations  (1.4)  simplify  as  follows: 

(k  -  1  )Nk(t)  =  tNk-y(t)  +  (k~  t)Nk .  ,(/  -  1).  (2.4) 

The  knot  sequence  t  =  IB  contains  just  two  points,  viz.,  the  points  0  and  1.  but  each 
with  infinite  multiplicity.  The  only  nontrivial  B-splines  for  this  sequence  are  those  which 
have  both  0  and  1  as  knots,  i.e.,  those  Blk  for  which  t,  =  0  and  =  1-  There  seems  to 
be  no  natural  way  to  index  the  entries  in  the  sequence  IB.  Instead,  it  is  customary  to  index 
the  corresponding  B-splines  by  the  multiplicities  of  their  two  distinct  knots.  Precisely. 

Bm  :=  £(•;  ,  1, . . . ,  1  ).  (2.5) 

1  times  v-ri  times 

With  this,  the  recurrence  relations  (1.4)  simplify  as  follows: 

B(p.„)(0  =  tB(M.„-n(i)  +  (1  -  (2.6) 

This  gives  the  formula 

£(„,„)(*)  =  (M  ^  ‘'Vi  -  for  0  <  t  <  1  (2.7) 

for  the  one  nontrivial  polynomial  piece  of  j?(fi  as  one  verifies  by  induction.  The  formula 
enables  us  to  determine  the  smoothness  of  the  B-splines  in  this  simple  case:  Since 
vanishes  identically  outside  0.1.,  it  has  exactly  v  -  1  continuous  derivatives  at  0  and 
p  -  1  continuous  derivatives  at  1.  This  amounts  to  v  smoothness  conditions  at  0  and 
p  smoothness  conditions  at  1.  Since  the  order  of  B(u^\  is  p  —  u  1.  this  is  a  simple 
illustration  of  the  generally  valid  formula 

^smoothness  conditions  at  knot  —  multiplicity  of  knot  -  order.  (2.8) 

For  fixed  p  -  u.  the  polynomials  in  (2.7)  form  the  so-called  Bernstein  basis  (for 
polynomials  of  degree  <  p  -  o)  and.  correspondingly,  the  representation 

P  =  YL  (-.9) 

H  —  i/—  h 

is  the  Bernstein  form  for  the  polynomial  p  f  r h .  In  CACJI).  it  is  more  customary  to 
refer  to  (2.9)  as  the  Bezier  form  (for  the  polynomial  p)  or  as  the  Bezier  polynomial  or 
even  the  Bernstein-Bezier  polynomial.  It  may  be  simpler  to  use  the  short  term  B-forin 
instead. 
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Figure  2.2  Bernstein  basis  of  degree  4 


3.  A  simplifying  assumption 

In  the  next  sections,  we  develop  the  basic  B-spline  theory  by  studying  the  spline  space 
Sk. t,  i.e..  the  collection  of  all  functions  s  of  the  form 


s 


B'ka' 


(3.1) 


for  a  suitable  coefficient  vector  a  -  (a,). 

In  practice,  the  knot  sequence  t  is  always  finite,  hence  so  is  the  sum  in  (3.1).  This 
often  requires  one  to  pay  special  attention  to  the  limits  of  that  summation.  Since  I  find 
that  distracting,  I  will  assume  from  now  on  that  the  knot  sequence  t  is  bi-infinite.  This 
can  always  be  achieved  simply  by  continuing  the  sequence  indefinitely  in  both  directions 
(taking  care  to  maintain  its  monotonicity)  and  choosing  the  additional  B-spline  coefficients 
to  be  zero. 

More  than  that,  I  will  assume  that 


/.rot  :=  lira  t,  -  :  oc. 

i  —  cx 


(32) 


This  assumption  is  convenient  since  it  ensures  that  every  r 
B-spline. 


IR  is  in  the  support  of  some 
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At  times,  it  will  be  convenient  to  assume  that 

t .  <  li+k  Vi  (3.3) 

which  can  always  be  achieved  by  removing  from  t  its  t  —  th  entry  while  i,  =  tt  +  k-  This 
does  not  change  the  space  Sk,t  since  the  only  Ar-order  B-splines  removed  thereby  are  zero 
anyway.  In  fact,  another  way  to  state  the  condition  (3.3)  is: 

Bik  #  0  Vt.  (3.3') 


4.  The  polynomials  in  t 

We  show  in  this  section  that  contains 

7r</t  :=  the  collection  of  all  polynomials  of  degree  <  k , 
and  give  a  formula  for  the  B-spline  coefficients  of  p  £  5r<fc. 

We  begin  with  Marsden’s  Identity: 


Theorem  4.  For  any  r  £  IR, 


(•  ~r)k-'  =]TBm(t), 

t 

(4.1a) 

with 

V'jfc(j')  :—  (<,.+  j  —  r)  •  •  •  (t i-rk-i  —  r)- 

(4.16) 

Proof 

sequence  a. 

W7e  deduce  from  the  recurrence  relation  (1.4)  that,  for  an  arbitrary  coefficient 

Z^BikCi  =  ^  5i,*-i((l  -  u ',fc)a,_,  *  u.tfcat). 

(4.2) 

On  the  other  hand,  for  the  special  sequence 

o.  :=  i\k{r)  :=  (t,-L i  -  r)  •••(*, +k-i  ~  t) 

(with  t  £  DR) 

.  we  find  for  Bl<k-\  *  0.  i.e.,  for  it  <  j  that 

(1  "  1 

*  kQ’  =  (f1  ~  ,«*'i>)(<i  -  t)  t-  Uik  ■  (ti+k- 1  -  r))v,,fc.-  |(r) 

=  (•  -  T)*i,k-l{T) 

(4.3) 

since  (1  -  ,)  is  the  straight  line  which  agrees  with  /  at  t,  and 

hence  must  equal  /  if,  as  in  our  case,  /  is  linear.  This  shows  that 

y]  BtkPikir)  ---  (■  -  t)  Bi,k-iWi,k-i(T). 


■r.  W/ .  y.  ■ 


10 


Carl  de  Boor 


hence,  by  induction,  that 

=  (•  -  r)*'1  ]TB.iV>ii(T)  =  (•  -  r)k~K 

since  v,i(r)  =  1  and  =  1  (see  (1.2)).  jji 

Remark  There  may  be  some  doubt  as  to  why  xl>x j  should  be  identically  equal  to  1. 
From  the  definition  (4.1b),  it  would  appear  that  0,j  is  the  product  of  no  factors,  hence, 
by  a  standard  agreement  concerning  the  empty  product,  equal  to  1.  This  is  the  definition 
appropriate  for  use  in  induction  arguments.  Indeed,  if  you  consider  the  coefficients  in  (4.2) 
for  k  —  2  directly,  you  get 

(1  -  W, -2)0i- 1.2(f)  +  u;,20i2(f)  =  (1  -  ^t2)  •  (<.  -  r)  +  •  (t,+  l  -  t) 

=  (•  -  0. 

which  agrees  with  (4.3)  for  this  case  if  we  set  ^ti(r)  =  1. 

Since  7  in  (4.1)  is  arbitrary,  it  follows  that  Sk,t  contains  all  polynomials  of  degree  <  k. 
More  than  that,  we  can  even  give  an  explicit  expression  for  the  required  coefficients,  as 
follows. 


By  dividing  (4.1a)  by  ( k  -  1)!  and  then  differentiating  it  with  respect  to  r,  we  obtain 
the  identities 


(It  -  1/)!  ^  ‘4  (k  -  I)! 


v  >  0, 


(4.4) 


with  Df  the  derivative  of  the  function  /.  On  using  this  identity  in  the  Taylor  formula 


(*-")! 


Dk~ 


W. 


valid  for  any  p  €  n<k,  we  conclude  that  any  such  polynomial  can  be  written  in  the  form 


P  = 


y ,  Bik^ikP  ■ 


(4.5a) 


with  \tk  given  by  the  rule 


k 

A,„/  :=  V 


l/-  1 


(-£>)"- V,fc(r) 

(A--  1 )! 


Dk-Vf{r). 


(4.5  b) 


Here  are  two  special  cases  of  particular  interest.  For  p  -  1.  we  get 


1  - 


(4.6) 
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since  Dk  V,*  =  (-1)*  *(Jk  —  1)!,  and  this  shows  that  the  B,k  form  a  partition  of  unity. 
Further,  since  Dk~2rlixk  is  a  linear  polynomial  which  vanishes  at 

*t  (*»+i  +  ’ '  •  +  ti+k-i) (k  -  1),  (4.7) 

we  get  the  important  identity 

p=^Bikp(rt)  Vp€ffi.  (4.8) 

t 

Remark  In  the  cardinal  case, 

V>«fc(r)/(*  -  1)!  =  ^  'fc**  Jj, 

while  in  the  Bernstein-Bczier  case, 

<t^M)  =  (-'•)“(!  -  0'  =  <-)'iw("  * ")• 


5.  The  pp  functions  contained  in  5jt,t 

In  this  section,  we  show  that  the  spline  space  Sk.t  coincides  with  a  certain  space  of 
pp(  :=  piecewise  polynomial)  functions. 

Each  s  t  Sk.t  is  PP  of  degree  <  k  with  breakpoint  sequence  t  since  each  Bxk  is  pp  of 
degree  <  k  and  has  breakpoints  f In  symbols. 

Sk.t  Q  *<k,t-  (5.1) 

But  Sk,t  is  usually  a  proper  subspace  of  n<k,t  since,  depending  on  the  knot  multiplicities 

(5.2) 

the  splines  in  Sk.t  are  more  or  less  smooth,  while  the  typical  element  of  n<k.t  bas  jump 
discontinuities  at  every  tx. 

For  the  precise  description,  given  in  Theorem  5  below,  of  the  smoothness  conditions 
satisfied  by  the  elements  of  Sk.t-  we  make  use  of  Marsden’s  Identity,  (4.1),  since  it  provides 
us  with  the  B-spline  coefficients  of  various  pp  functions  in  Sk.t,  as  follows. 


(4.16)2- 


(4.16)b 


Carl  de  Boor 


i 

V 

\ 

^  i  —  1  ^  t  -f  2  ^ t  +  3 \  ^i+4 


Figure  5.1  Bl4  and  ipi4\  note  the  double  knot 

Since  Btk(tj)  ¥■  0  implies  rl-\k{tj)  =  0  (see  Fig.  5.1),  the  choice  r  =  tj  in  (4.1)  leaves  only 
terms  with  support  either  entirely  to  the  left  or  else  entirely  to  the  right  of  t}  \  see  Fig.  5.2. 
This  implies  that 

*■>>■ 

with 

Qj.  :=  max{o,0}  (5.4) 

the  positive  part  of  the  number  a.  More  than  that ,  since  B,k  (f  j )  t  0  implies  Dv~  Vifc(fy)  = 
0  in  case  v  <  £ty  the  same  observation  applied  to  (4.4)  shows  that 

(•-*,)$""  €  Sfc.t  for  1  <  */<  #<j.  (5.5) 


Theorem  5.  The  space  S*it  coincides  with  the  space  5  of  all  piecewise  polynomials 
of  degree  <  k  with  breakpoints  t,  which  are  k  1  -  #/,  times  continuously  differentiable 
at 

Proof  Assume  without  loss  of  generality  (see  Sec.  3)  that 

<  t i  +  fc  '*!• 

It  is  sufficient  to  prove  that,  for  any  finite  interval  I  j a.b.,  the  restrict  ion  5  /  of  the  space 
5  to  the  interval  /  coincides  with  the  restriction  of  S*it  to  that  interval.  The  latter  space 
is  spanned  by  all  the  B-splines  having  some  support  in  I.  i.e.,  all  f?,*  with  (/,.  *1 r  /  *  0. 

The  space  Sj  has  a  basis  consisting  of  the  functions 

(■  -  a)  ,  v  =  (■  -  v  1.  ■  •• .  for  a  <  l,  <  b.  (5.6) 


Figure  5.4  (a)  The  six  quadratic  B-splines  for  the  case  of  one  simple  and 
one  double  interior  knot:  and 
(b)  the  corresponding  truncated  power  basis. 

On  the  other  hand,  the  dimension  of  S‘j,  i.e.,  the  number  of  functions  in  (5.6),  equals  the 
number  of  B-splines  with  some  support  in  1  (since  it  equals  k  +  $Z0<t,<f>  hence  is 
upper  bound  on  the  dimension  of  (S*,t)|/-  This  implies  that  equality  must  hold  in  (5.7). 
which  is  what  we  set  out  to  prove,  jl! 

Remark  The  argument  from  Linear  Algebra  used  here  is  the  following:  Suppose 
that  we  know  a  basis,  (/i,  /2>  •  •  • ,  fn)  say,  for  the  linear  subspace  F ,  and  that  we  further 
know  a  sequence  (g\,g2i  •  •  •  •  9m)  whose  span,  G  say,  contains  each  of  the  /,.  Then,  of 
course,  F  C  G  and  so 

n  —  dim  F  <  dim  G  <  m. 

If  we  now  know,  in  addition,  that  n  —  m,  then  necessarily  F  —  G.  Moreover,  then 
riecesarily  dim  G  ~  m.  i.e..  the  sequence  (<7i-?2*  •  •  •  i9m)  must  be  linearly  independent 
(since  it  then  is  minimally  spanning  for  G).  In  our  particular  situation,  this  last  observation 
implies  that  the  set  of  B-splines  having  some  support  in  /  must  be  linearly  independent 
over  /.  We  pick  up  on  this  in  the  next  section. 
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6.  ‘B’  stands  for  ‘BASIC' 

In  this  section,  we  discuss  the  basis  property  of  the  B-splines.  as  a  consequence  of 
Theorem  5  and  its  proof. 

From  the  Remark  following  Theorem  5,  we  obtain  the  following  sharpening  of  Theorem 
5. 


Theorem  6.  Let  I  [a,  b[  be  a  finite  interval.  Then  the  restrictions 


{Btk\j  ■  Bik\i  /  °}  (6-1) 

of  those  B-splines  which  have  some  support  on  I  form  a  basis  for  the  space  of  pp  functions 
of  degiee  <  k  on  I  with  breakpoints  {tt  :  a  <  t,  <  b}  and  which  are  k-  1  -  continuously 
differentiable  at  each  of  their  breakpoints  t,. 

We  conclude  that  the  number  of  smoothness  conditions  at  a  knot  t,  guaranteed  to  be 
satisfied  by  every  spline  in  Sk,t  equals  k  -  #tt.  This  proves  the  formula 

^smoothness  conditions  at  knot  +  multiplicity  of  knot  —  order  (2.8) 

cited  earlier  (in  connection  with  the  Bernstein-Bezier  form). 

It  is  worthwhile  to  think  about  this  the  other  way  around.  Suppose  we  start  off  with 
a  partition 

a  =:  <  s2  <*■•<&<&-*- 1  :=r  * 


of  the  interval  1  a.  b  and  wish  to  consider  the  space 


*<k,‘ 

of  all  pp  functions  of  degree  ;  k  on  /  with  breakpoints  £,  which  satisfy  ux  smoothness 
conditions  at  £,.  i.e.,  are  u,  1  times  continuously  differentiable  at  £,,V/.  Then  a  B-spline 
basis  for  this  space  is  provided  by  (6.1).  with  the  knot  sequence  t  constructed  from  the 
breakpoint  sequence  £  in  the  following  way:  To  the  sequence 


L) 


U'i-  •  •  ■  •  sC2;  ,  s3 - ,€t 

i/_-  times  i,-  times  v,  times 


(6.2) 


adjoin  at  the  beginning  k  points  <  a  and  at  the  end  k  points  >  6.  While  the  knots  in 
(G.2)  have  to  be  exactly  as  shown  to  achieve  the  specified  smoothness  at  the  specified 
breakpoints,  the  2k  additional  knots  are  quite  arbitrary.  They  are  often  chosen  to  equal  u 
resp.  6.  and  this  has  certain  advantages  (among  other  things  that  of  simplicity).  With  such 
a  choice,  it  is  necessary  to  modify  the  definition  (1.1)  so  as  to  include  the  right  endpoint. 
6.  into  the  support  of  the  rightmost  nontrivial  In  other  words,  if  n  is  such  that 
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then 

Bni(t)  :=  Xn(t) 

This  ensures  that,  in  evaluating  a  spline  or  its  derivatives  at  6,  we  obtain  the  limit  from 
the  left. 

The  identification  of  S*.t  with  a  certain  space  of  pp  functions  allows  the  following 
conclusions  of  importance  in  calculations  to  be  discussed  later. 

Corollary  1.  If  <t  <  <,  +  *_!.  then  the  derivative  of  a  spline  in  Sk>t  is  a  spline  of 
degree  <  k  -  1  with  respect  to  the  same  knot  sequence,  i.e.,  DSk< t  £ 

Proof  By  assumption,  #t,  <  k ,  hence  the  pp  functions  in  Sk<t  are  continuous, 
therefore  differentiable  (if  we  accept  a  possible  jump  at  f,  in  the  derivative  Ds  of  s  6  §k,t 
in  case  -  k  -  1).  Further,  such  a  derivative  Ds  is  pp  of  degree  <  k  -  1  and  satisfies 
k  ~  Mtt  -  1  smoothness  conditions  at  tt ,  hence  belongs  to  Sk~  l.t,  by  Theorem  5  or  6.  ||| 

Corollary  2.  If  t  is  a  refinement  of  the  knot  sequence  t,  then  Skt t  C  Sk 

Proof  Since  t  is  a  refinement  of  t,  i.e.,  contains  entries  in  addition  to  those  of  t, 
the  pp  functions  in  Sk<t  satisfy  ail  the  conditions  which,  by  Theorem  5  or  6,  characterize 
the  pp  functions  in  Sk  £.  (But  the  converse  does  not  hold,  since  the  pp  functions  in  Sk  £ 
may  have  more  breakpoints  and/or  may  be  less  smooth  at  some  breakpoints  than  the  pp 
functions  in  S*,t.)  I’- 

These  corollaries  point  out  that  it  should  be  possible,  in  principle,  to  compute  from 
the  B-spline  coefficients  of  a  spline  in  Skj  the  B-spline  coefficients  of  its  derivative  and  its 
B-spline  coefficients  with  respect  to  a  refined  knot  sequence.  To  carry  out  such  calculations, 
though,  we  need  a  means  of  expressing  the  B-spline  coefficients  of  a  spline  in  terms  of  other 
information,  such  as  its  values  and  derivatives  at  certain  points.  If  the  spline  happens  to  be 
a  polynomial,  then  such  a  formula  is  provided  by  (4.5).  We  show  in  the  next  section  that 
the  same  formula  works  for  any  spline  (provided  we  are  willing  to  restrict  the  parameter 
t  suitably). 


1,  if  ti<t<b 
0,  otherwise. 


(6.3) 


7.  The  dual  functionals 


In  this  section,  we  prove  that  the  formula  (4.5)  for  the  B-spline  coefficients  of  a 
polynomial  is  valid  for  an  arbitrary  spline  provided  we  restrict  the  parameter  t  in  the 
definition 

L  -  (FTT)— 


(4.56) 


to  the  support  of  Blk. 
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For  this,  we  agree,  consistent  with  (1.1b),  that  all  derivatives  in  (4.5b)  are  to  be  taken 
as  limits  from  the  right  in  case  r  coincides  with  a  knot  (except,  perhaps,  when  r  is  the 
right  endpoint  of  the  interval  of  interest,  see  (6.3)). 


Figure  7.1  The  three  polynomials,  p/_  2,P:-i-Pi,  which  agree  with  some 
quadratic  B-spline  on  the  knot  interval 


Theorem  7.  If  r  in  definition  (4.5b)  of  A,*  is  chosen  in  the  interval  jt,-, then 


^  B j fc Q j  )  —  a,- 


Proof  YVo  prove  that,  under  the  given  restriction. 


^ikBjk  —  — 


1.  if  J  -  J\ 

0,  otherwise. 


Assume  that  ?  t  ,//,  j;C  r*  .  Then  (7.2)  requires  proof  only  for  j  =  l  -  k  t  1, - 

since,  for  ail  other  j.  t  ^  j  and  B-,k  vanishes  identically  on  .  hence  also  \tkB}k  -  0 

for  each  of  the  remaining  j' s.  let  p,  be  the  polynomial  which  agrees  with  B, k  on 
Then 

^ikBjk  ~  ^rkPj  ■ 


On  the  other  hand, 


P :  ~  ^  P i  ^ikPj 


i  =  l-  k- I 


since  this  holds  by  (4.5a)  on  [</, f/+i[.  This  forces  A tkp:,  hence  \xkB}k,  to  equal  6,;  for 
i,j  =  /  -  k  +  1, . . . ,/,  since,  by  Theorem  6  or  directly  from  the  fact  that  (4.5a)  holds  for 
every  p  €  n <*,  the  sequence 

Pl-k-^u"  ■  iPl  (7-4) 

is  linearly  independent.  j| 1 

Remark  The  argument  used  here  is  that,  for  a  linearly  independent  sequence 
(/i<  •  •  • ,  /*).  the  only  way  the  equation 

n 

ft  -  ^  / )atj 

can  hold  is  for  atJ  to  equal  1  for  t  =  j  and  zero  otherwise.  Further,  the  linear  independence 
of  the  sequence  (7.4)  follows  from  the  validity  of  (4.5a)  for  every  p  €  KKk  since  that  implies 
that  the  ^-sequence  (7.4)  is  spanning  for  the  ^-dimensional  space  it<k-  It  also  follows 
from  Theorem  6  with  /  =  ti,t[+ il. 

The  two  sequences,  (B,*)  and  (A^*),  are  said  to  be  bi-orthonormal  or  dual  to  each 
other  because  they  satisfy  (7.2).  For  this  reason,  the  linear  functionals  A,*  are  at  times 
referred  to  as  the  dual  functionals  for  the  B-splines. 

We  exploit  the  simple  formula  (7.1)  for  the  a'-th  B-spline  coefficient  of  a  spline  in 
subsequent  sections,  in  order  to  derive  algorithms  for  differentiation  and  knot  insertion  and, 
ultimately,  to  derive  statements  about  the  condition  and  the  shape-preserving  property  of 
B-splines. 


100  101 


Figure  8.1  The  power  coefficients  of  these  two  very  different  linear  poly¬ 
nomials  differ  by  only  0.1 '7. 
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8.  Condition 

The  condition  of  a  basis  measures  how  closely  relative  changes  in  the  coefficients  are 
matched  by  the  resulting  relative  changes  in  the  element  represented.  The  closer  the  match, 
the  better  conditioned  the  basis  is  said  to  be.  For  example,  the  power  basis  1  .  is 

not  a  good  way  to  represent  polynomials  if  we  are  interested  in  a  positive  interval  [a,  A- 
with  a/fe  close  to  1.  If,  e.g.,  'a,  b  =  [100. 101  .  then  a  0.1  %  change  in  the  power  coefficients 
of  the  straight  line  p  :  t  — >  t  -  100  can  change  its  behavior  on  '  1 00, 101  i  by  100%;  see 
Fig.  8.1. 

If  we  use  the  appropriately  shifted  power  basis,  e.g.,  write  p  in  the  form  p(t)  — 
a  -  j3(t  -  100).  then  a  .1%  change  in  the  coefficients  a, 0  of  this  form  produces  a  .1% 
change  in  the  polynomial  on  the  interval  [  100, 101  i .  The  appropriately  shifted  power  basis 
is  often  much  better  conditioned  than  the  power  basis.  In  this  section,  we  discuss  briefly 
the  condition  of  the  B-r  pline  basis. 

This  requires  us  to  bound  the  spline  in  terms  of  its  B-spline  coefficients  and  the  B- 
spline  coefficients  in  terms  of  the  spline.  The  first  turns  out  to  be  easy,  while  the  second 
requires  some  work.  Precisely,  we  are  looking  for  constants  m  >  0  and  M  for  which  the 
inequalities 

m  max  !a,i  <  max  ;  V'  B,*(f)at  <  .W  maxjO,  (8.1) 

i  t  i 

l 

hold  regardless  of  what  the  coefficient  vector  a  =  (a, )  might  be.  Since  the  B-splines  are 
nonnegative  and  sum  to  1  at  any  point,  we  have 

'X>  k{t)at  <  Y'  Bik{t)\a,  <  Y  Bt*(f)  max  ,a,:  =  max  a,  , 


hence  the  second  inequality  always  holds  with  M  =  1.  For  the  first  inequality,  we  have  to 
work  a  little  harder. 

Set  s  Bikai-  We  know  from  Theorem  7  that 


a ,  =  ~ 


{-D)u  V’.feN 
'  (*-!)!  °  "M 


with  t  some  point  which  we  can  freely  choose  in  the  interval  7,. .  We  now  bound  this 
sum  in  terms  of  max*  's(t)'.. 

Suppose  that  r  €  i:C  (,,<,+^1.  Then,  for  some  const^  depending  only  on  k , 

and  for  all  p  C  n <*.  and  all  j. 

D: p(t)  <  const*.  (A?*)  max  p(t)  .  (8.3) 

The  existence  of  such  a  const*  follows  for  the  case  A t\  =  1  froin  the  fact  that  7r<*  is 

finite-dimensional,  and  from  this  it  follows  for  arbitrary  At,  by  scaling.  Since  ,s  agrees 

with  some  polynomial  of  degree  <  k  on  .  we  conclude  that 

Djx(t)  <  const*  ( A//)  ’ ;  max  *(f)  .  (8.4) 

t, 
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On  the  other  hand,  t>,*  =  (tx+\  (tx^k~  \  ~  )  is  also  a  polynomial  of  degree  <  k,  and 

max  v, k(f)  <  const i  A//-  k~  1  (8.5) 

U<t  S'i-i 

lor  some  ronst^  which  depends  only  on  k  and  with  //-,</-  +  i;  a  largest  interval  of  that 
form  in  i/,, .  Therefore  we  choose  1  =  1’  and  then  obtain,  from  (8.3)  with  p  =  tjjlk  and 
from  (8.4),  the  bound 

Dl'~lxi.'ik(r)Dk-i'S{r)\  <  (const*)2const'fc  max  |s(t);. 

t, 

Now  sum  these  bounds  over  v  and  divide  by  (k  -  1)!  to  obtain 

flt!  ~  <  const  max  «(/)•', 

with  const  depending  only  on  k. 

We  have  proved  the  following 

Theorem  8.  There  exists  a  constant  Dk  depending  only  on  k  so  that,  for  all  knot 
sequences  t  and  all  .<s  £  Sk  t,  and  for  all  t, 

Xlks\  <  Dk  max  |,«(/)i.  (8.6) 

«.<<<«.** 

The  best  value  for  Dk  is  not  known  exactly  but  there  is  strong  numerical  evidence 
that  Dk  ~  2*_1.  If  we  only  consider  cardinal  splines,  i.e..  only  uniform  knot  sequences, 
then  the  best  value  for  Dk  is  known  to  be  less  than  (tt/2)* 

Corollary.  The  inequalities  (8.1)  hold  with  m  =  1  Dk  and  M  =  1. 

9.  Evaluation 

In  this  section,  we  discuss  the  use  of  the  recurrence  relations  (1.4)  for  the  evaluation 
of  a  spline 

s=^B,kax  (9.1) 

I 

from  its  B-spline  coefficients  (a,). 

We  already  observed  in  (4.2)  that  the  recurrence  relations  imply 

g  =  'S  _  B,ka,  -  ^  Buk  .  i a,1  .  (9.2) 

/  I 

with 

n,1  :  -  (1  a,  |  -  ~,ka,.  (9.3) 
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Note  that  a  *  is  not  a  constant,  but  is  the  straight  line  through  the  points  (t,,a,_|)  and 
(f i  -  1. a,).  In  particular,  fl,1  (t)  is  a  convex  combination  of  a,-  j  and  a,  ift,  <  t  <  j. 

After  k  ~  1-fold  iteration  of  this  procedure,  we  arrive  at  the  formula 


which  shows  that 


k- 1;  ,  .  i 

•s  --  a ,  on  ir 


Algorithm  9.  From  given  constant  polynomials  a(  '  o,,  i  =  j  -  k -t  l,. . . ,  j, 

(which  determine  .«  :=  on  +  generate  polynomials  a\r  ,r  =  -  1. 


by  the  recurrence 


(1  “  ~Yk- r)n|r. ,  -  rfl/1.  j-*:-r+l<i<  j.  (9.4) 

Then  s  -  on  ljjj.,1'.  Moreover,  for  /,  <  £  <  £J  +  1,  the  weight  in  (9.4) 

lies  between  0  and  1.  Hence  the  computation  of  s(£)  =  af  1  (t)  via  (9.4)  consists  of  the 
repeated  formation  of  convex  combinations. 

In  the  cardinal  case  (see  Sec.  2,  esp.  (2.2-4)),  the  algorithm  simplifies,  as  follows.  Now 
•<  ==  V  A*(-  -  i)ax  =  -  *>!’  /(*  -  !)• 


Hence 


Q,  ;  :=(»  +  fc  --  1  -  -)at-i  +  (•  -  !>>• 


(fc  -  1)!  on  j.  j  4  1  , 


«r  :={;~k-r- 1  m-  *k  >-* 


r  *-  1  ./• 


In  the  Bernstein-Bezier  case  (see  Sec.  2.  esp.  (2.5-9)).  all  the  nontrivial  weight  func¬ 
tions  -c,.a  r  are  the  same.  i.e.. 

r(0  -  <• 


Thus,  for 


we  get 


.  n.u.,.  • 


£4  *  1/  -  /l 


-  -  -  .  A.s-y.l  -V-V* V-"X  \ -V  -\  0 

•’**'**  a  JSEt*  -*  **  \  *.  « %*  •  * *  ■ « ,  ■*  „  ^  j-  ,*  *.*  *  *  ■ 

--  1  « -mV.  j- V V'\/ .• 


/  V  V  *„• 


•  -  .  •  .A  A  ^  •*,  .*  , 

*,»  Vs/  >  •  .  •  .  * 

*  -VvV 
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■S  “  a(0,0)  on  |0,  1  , 

with  (9.4)  ib 

=  (1  -  t)a(n +  i,„)  +  *0(„,„-ir  H  +  ^  -  r:  r  =  /t  -  1 . 0. 

This  is  de  Casteljau’s  Algorithm  for  the  evaluation  of  the  B-form. 


10.  Differentiation 

In  this  section,  we  derive  a  formula  for  the  B-spline  coefficients  of  the  derivative  of  a 
spline  in  terms  of  the  B-spline  coefficients  of  the  spline. 

By  Corollary  1  to  Theorem  6,  the  derivative  Ds  of  a  spline  s  €  Sk,t  is  again  a  spline 
with  the  same  knot  sequence  but  of  one  order  lower.  This  means  that,  by  Theorem  7,  we 
can  compute  its  B-spline  coefficients  ( a[ )  by  the  formula 

a[  —  ^i.k- 1  {Ds) 


provided  we  use  r  6  jt, j!. 

To  relate  o'  to  a,  we  express  X^k-xD  as  a  linear  combination  of  the  functionals  A,*, 
making  use  of  the  fact  that  A,*,  depends  linearly  on  rpik~  recall  the  definition 


(  D)"  ‘V,»(0  Dk-»f{Ty 


V  -  1 


and  that 

These  facts  imply  that 


(*-!)! 

(f  i  -  k-  1  ~  f  i )  t't.A:  -  t  “  V\k  ~  t-'t  -  1  ,k- 


(4.56) 


(10.1) 


( A,<  A i  - 1  .it )/(r ) 


i/..  i 


(*-  1  )• 


<r-  1 


--  i 


(A  1 )! 


the  last  equality  by  (10.1)  and  since  Dk  \  0.  On  the  other  hand,  directly  from 

t  he  definition  (4.5b), 
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Comparison  of  these  two  displays  shows  that 

k  -  1 


A, 


,k-\D  =  - - —  (A,fc  -  A,_lt*  ). 

li-t  k~l  -  ' 


(10.2) 


Assuming  that  Bt,k-\  0,  i.e.,  that  t,  <  tx+k-i.  we  can  choose  r  c  (tx,tx  +  k- i)  - 
(/,_!, n  (l,,t,^k).  This  yields 

Algorithm  10.  Compute  the  coefficients  for  Y^a[Bx,k- 1  :=  D^axBxk  by 

ai_a,~1  :r  ■*  '  *  (10.3) 


a,  - 


'  (f,+fe_i  -  «*)/(*  -  1) 


,  if  1 1  <  1  —  l  • 


Remark  What  happens  when  t,  =  In  this  case,  BXik-i  -  0,  hence  there 

is  no  need  to  calculate  ax.  To  be  precise,  in  this  case,  the  spline  s  =  Y^,  Blkai  may  not 
even  be  continuous  at  tx.  therefore  ( Ds)(tx )  makes  no  sense.  On  the  other  hand,  the 
left  and  the  right  limit,  (Ds)(tx~)  and  (Ds)(t{+),  always  make  sense,  and  the  algorithm 
would  provide  all  the  a'^s  needed  for  their  calculation.  By  applying  the  algorithm  to  the 
particular  coefficient  sequence  a  -  (6X]).  we  obtain  the  formula 


DB,k  = 


k  -  1 


1 1  ^  k  —  l  f  t 


Bi,k- 1  - 


f  t+/fc  C-*-i 


B 


i+l,k-l- 


(10.4) 


In  terms  of  the  alternative  notations  (1.9)  for  B-splines,  this  reads 


DNik  =  Mi,k-i  ~  Mt-t-i.fc-i- 


Since  Nxk  -  1,  this  implies  that 


J  Mx,k-i  =  1  (10.5) 

and  so  indicates  why  the  particular  normalization 

A iik  - - -S»fc 

U~k  -  U 

is  of  interest. 

In  the  cardinal  case,  (10.3)  reduces  to 

a[  =  a,  -  a,_  ,  Va,  (10.3) zz 

and  (10.4)  reads 

DNk  =  A\_i  -  .Vfc-,(-  -  1)  (10.4)2 
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On  integrating  this  formula,  we  obtain 

Nk{t)=  f  Nk-i(r)dT  (10.6) 

Jt-  1 

since  both  sides  of  (10.6)  vanish  for  negative  t.  In  terms  of  the  convolution  product 

(/*<?)( 0  :=  J  f(t  ~  T)9(r)dT 

of  two  functions  /  and  g,  this  gives  the  important  formula 

tf*  =  N*-i*.Vi(-  +  l).  (10.7) 

This  shows  that  Nk{-  +  k  -  1)  is  the  it- fold  convolution  product  of  N j,  i.e., 

Nk(-  +  k  -  1)  =  AT)  *  Ni  *  . . .  *  JVj . 


k  terms 


In  the  Bernstein- Bezier  case,  we  get 


with 


p+v—h  n+v=h-l 


a{ti,u)  —  (ll  +  V  +  1)  l,i/) 


(10.3) 


IB 


11.  Knot  insertion 

In  this  section,  we  discuss  the  most  important  CAGD  contribution  to  (univariate) 
spline  theory,  viz.,  the  idea  of  knot  insertion  (a.k.a.  subdivision).  Since  the  spline  order, 
k ,  will  not  change  in  this  section,  we  will  usually  suppress  it  and  write  Bt  instead  of  /?,•*, 
ipt  instead  of  etc. 

Simply  put,  knot  insertion  involves  rewriting  a  given  spline  as  a  spline  with  a  refined 
knot  sequence,  as  can  always  be  done  by  Corollary  2  of  Theorem  6.  Such  a  calculation  is 
worthwhile  since  the  B-spline  coefficients  are  nearly  equal  to  values  of  the  spline  at  known 
points,  and  this  is  more  nearly  so  when  the  knots  are  closer  together.  Here  is  the  precise 
statement. 

Theorem  11.  If  the  spline  s  -  Yi,  B,a ,  is  continuously  differentiable,  then 

la,  -  «(<*)  <  const  t  2  sup  D2s{t )  .  (11.1) 
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f (tt  +  i  +  *t+2  *  *«-+fc-i)/(^  _  1) 


It)  :=  sup(*t,!  -  <,). 


Proof  Recall  from  Sec.  8  that 


la,1  =  !A,sl  <  const  max  is(t)j. 


Further,  recall  from  Sec.  4  (esp.  (4.8))  that 

A iP=p(t't)  6  7T i . 

Thus,  choosing,  in  particular,  p  :=  s(t’)  +  (•  -  t[)Ds(t j),  the  linear  Taylor  polynomial  for 


s  at  /  we  get 


1 0,  -  «(/;)!  =  ia,  -  p(<;);  =  jA,(s  -  p)l  <  const  max  \(s  -  p)(t)! 


v  const  k — max  !£>2s(f)!. 

8  +  k 


Figure  11.1  A  cubic  spline  and  its  control  polygon.  The  end  knots  are 
quadruple. 


This  suggests  consideration  of  the  control  polygon  associated  with  the  representa¬ 
tion  B,a,  of  the  spline  s  as  an  element  of  •  This  control  polygon  will  be  denoted 
bv 

C„.«. 
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It  is  the  broken  line  or  piecewise  linear  function  with  vertices  Px  :=  For,  the 

theorem  implies  that  the  control  polygon  will  be  close  to  s  if  it:  is  small.  Here  is  the 
precise  statement. 


Corollary.  Let  C„. t  be  the  control  polygon  associated  with  the  representation 
!C«  ^>ai  °f  the  continuous  spline  s  as  an  element  of  St.  Then 


sup|s(f)  -  CQ,t(t);  <  const  t  2 sup  D2s{t)>. 


(11.2) 


Proof  Let  tx  <  t  <  <t  J_,  and  let  p  be  the  linear  polynomial  which  agrees  with  s  at 
t'x  and  /t  +  1.  Then 

i*(0  -  p(0!  <  -  tf'|2/8  max  I D2s(t)\, 


while 


P(0  S  rnax{is(t1 )  -  a,j,  -  at4.i ,}  <  const  it j 2  max  !L>2s(r)j 


by  the  theorem. 


figure  11.2  The  control  polygon  of  Fig.  11.1  and  three  midpoint  refine¬ 


ments. 


This  shows  that  the  control  polygon  C„.t  converges  to  the  spline  5  as  we  refine  the 
knot  sequence  t.  Since  the  typical  graphical  equipment  only  draws  broken  lines,  anyway, 
this  makes  it  attractive  to  construct  refined  control  polygons  for  a  spline. 

for  this,  we  need  to  know  how  to  compute,  from  its  B-spline  coefficients  o,  as  an 
element  of  St.  the  B-spline  coefficients  a,  for  the  spline  .s  with  respect  to  a  refined  knot 


mmm 
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1 


sequence  t.  By  Theorem  7,  this  is  a  question  of  comparing  the  corresponding  A,  with  A,. 
Since  the  dual  functional 


i/=  i  ' 


(4.56) 


depends  linearly  on  this  requires  nothing  more  than  to  express 

b’t  =  (t t  +  k -  l  ~  ‘ ) 

as  a  linear  combination  of  the 

This  is  particularly  easy  when  t  is  obtained  from  t  by  adding  just  one  knot,  say  the 
point  i.  Then 


l  t  <  tt , 


hence  there  is  some  actual  computing  necessary  only  for  t,  <  t  <  1-  For  this  case, 

att’.-i  t-  Mi  -  (*,  +  i  -  •)•••(<, ^fc_2  ~  -  •)  +  &{ti+k~\  -  -)j 


provided  o(/,  -  •)  +  d(t,  -  •)  =  (t  -  ),  i.e., 

a  =  l-  u,,(f)  and  jH  =  tu,(t). 

Since  f,  --  /,  <  t  <  tt  +  k-\  ~  tx+k,  we  can  choose  r  in  the  definition  (4.5b)  in  the  interval 
(ttJirk)  -  U,  +  t)  r;  This  proves 

Algorithm  11.  If  the  knot  sequence  t  is  obtained  from  the  knot  sequence  t  by 
addition  of  the  point  t.  then  the  coefficients  a,  for  the  spline  s  with  respect  to  the  refined 
knot  sequence  are  given  by 


{®ii  ^  h  **  -]  ^  fi 

(1  —  w,(r))Gj_i  u.'I(f)flt,  iff,  <  t  <  f  i  j-fc-  i  J 

Oi  l,  if  t  <  t,. 


(11.3) 


Observe  that  u?i(i)  t  0.  ll  when  t,  <  i  <  t,+k-  i-  and  thus  the  coefficients  a  are  convex 
combinations  of  the  coefficients  a. 

Th  is  algorithm  has  the  following  very  pretty  graphical  interpretation. 

Corollary.  The  refined  control  polygon  C-  f  can  be  thought  of  as  having  been  ob¬ 
tained  by  interpolation  at  its  vertices  to  the  original  control  polygon  Ca-t,  i.e.. 


Ca.i(r>)  =  Ca.*(0  for  all/. 


(11.4) 


'SSXSm 
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Proof  Consider  the  straight  line  p:t^t.  It  is  a  spline  and,  by  (4.8), 

p  = 

t 


i.e.,  (<j )  is  its  B-spline  coefficient  sequence  with  respect  to  the  knot  sequence  t.  In  partic¬ 
ular,  it  is  its  own  control  polygon,  i.e.,  Ct-,t  =  P,  regardless  of  what  the  knot  sequence  t 
might  be.  This  implies  that  (11.3)  also  holds  with  every  a  replaced  by  r.  Pj 

This  says  that  the  point  P,  :=  (t}.aj)  lies  on  the  segment  and  cuts  this 

segment  in  the  ratio  (t  -  t3)  :  [t3+k- i  —  This  is  illustrated  in  Figure  11.3  for  the  control 
polygon  of  Figure  11.1. 


tj-l  t  tj  ^  l 

0—9 - 0 


Pj 


-i 


1  0 

i  o 

tj  —  2  () 

Figure  11.3 


Insertion  of  t  =  2  into  the  knot  sequence 
t  =  (0,0, 0,0, 1,3, 5, 5, 5, 5),  with  k  =  4. 


Q  tj 

<rt 


6  ^-^3 


1 

1 


If  r  :=  #f  <  k  —  1,  then,  after  just  (it  —  1  -  r)-fold  insertion  of  t.  we  obtain  a  knot 
sequence  t  in  which  the  number  i  occurs  exactly  k  -  1  times.  This  means  that  there  is 
exactly  one  B-spline  for  that  knot  sequence  which  is  not  zero  at  t.  Hence  it  must  equal  1 
at  t  and  its  coefficient  must  provide  the  value  of  s  at  t.  This  makes  it  less  surprising  that 
the  calculations  in  Algorithms  9  and  11  are  identical. 


Figure  11.4  The  cubic  spline  and  its  control  polygon  from  Figure  11.1  and 
the  sequence  of  control  polygons  generated  by  three-fold  inser¬ 
tion  of  the  same  knot.  (The  finest  control  polygon  differs  from 
its  predecessor  only  by  an  additional  vertex  point.) 


Figure  11.5  Conversion  to  B-net  by  (k  -  2)-fold  insertion  of  each  knot 

Conversion  to  B-net  Let  t'  be  the  refined  knot  sequence  which  contains  each  of 
the  knots  in  t  exactly  k  -  I  times.  Then  each  corresponding  B-spline  B' k  is  nonzero  on 
just  one  knot  interval,  hence  coincides  there  with  a  properly  shifted  and  scaled  element  of 
the  Bernstein  basis.  The  k  B-spline  coefficients  o'  associated  in  this  way  with  a  knot 
interval  therefore  provide  the  coefficients  in  the  B-form  for  the  polynomial  with  which  the 
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spline  agrees  on  that  knot  interval.  The  coefficient  sequence  (a[),  or  the  control  polygon 
Ca'  t<.  are  called  the  B-net  for  the  given  spline.  It  can  be  obtained  by  inserting  each  knot 
l,  of  the  spline  k  —  1  #<,•  times.  The  process  can  be  speeded  up  slightly  by  inserting 

first  every  other  knot,  and,  in  a  second  round,  inserting  the  remaining  knots.  The  latter 
insertion  process  is  then  entirely  local  and  depends  only  on  the  ratio  of  the  two  knot 
intervals  containing  the  knot  being  inserted. 

While  the  formulas  do  simplify  for  the  cardinal  case,  they  are  not  of  much  use  in  that 
form  since  insertion  of  one  knot  into  the  sequence  t  —  TL  would  destroy  the  uniformity  of 
the  knot  sequence.  But.  it  makes  good  sense  to  develop  formulas  for  inserting  the  same 
number  of  uniformly  spaced  knots  into  every  interval  jt,t  +  1  since  this  produces  again  a 
uniform  knot  sequence.  Because  of  its  practical  importance,  we  treat  this  case  separately, 
in  the  next  section. 


12.  Knot  insertion  for  cardinal  splines 

In  this  section,  we  consider  knot  refinement  for  cardinal  splines,  i.e.,  splines  with 
a  uniform  knot  sequence.  Here  it  is  desirable  to  have  the  refined  knot  sequence  again 
uniform.  We  restrict  attention  to  the  case  that  the  given  knot  sequence  is  t  —  TL.  This  is 
no  real  restriction  since  an  arbitrary  uniform  knot  sequence  can  always  be  written  in  the 
form  a  -  3TL  for  appropriate  scalars  a  and  3.  and  if  s  is  a  spline  with  that  knot  sequence, 
then  s(o  --T-  3  )  is  a  spline  with  the  knot  sequence  TL. 

If  we  insert  rn  -  1  uniformly  spaced  knots  into  every  knot  interval  of  K,  then  the 
refined  knot  sequence  is  t  -  TL  m.  The  corresponding  B-splines  B,  are 

Bx  =  .V *(•  -  t), 

with 

AT(0  N * (mt) 

an  appropriately  scaled  version  of  the  standard  cardinal  B-spline  A*..  This  makes  it  trivial 
to  determine  a,  in  case  k  —  1.  Since 

V,  -  A',  A’,  (■  -  rn-]),  (12.1) 

we  fmd  for  this  case  that 


ar ,+j  -  a,  for  j  =  0. 


rr,  -  1. 


The  formula  for  general  order  k  is  obtained  from  this  with  the  aid  of  the  convolution 
formula 


AT  -  A 


k~  -'k  1 


.V, 


1 


(10.7) 


Irorn  Sec.  10.  as  follows.  We  define 


t 


IT* 


r  1 


k. 


(12.2) 
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Then  a,  =  a,*,  and,  from  (10.7)  and  (12.1), 

m 

«r+ 1  =  Sr  *  iVj(-  +  1)  =  (  ^(*  ~  *>»«•)  *  £1  '  j) 

*  >=I 

m 

-EE  *,(■-. )*.V|(--  »  a., 

^  _  j  *  . .  1  “s/*  1  —  —  " 

.Vr+j(-  -  i  -  j  +  l)/m 

m 

=  5H  ^*r+l(*  —  *)/m  ,r. 

«  J=1 

Here,  we  have  used  the  following  consequence  of  the  convolution  formula  (10.7): 

•Vr_i(-  -  a)  *  Arj(-  -  0)  =  J  Nr-i{m(-  -  t)  -  a) A ^(mr  -  3 )dr 

=  J  JVr_i(m  ■  -<?  +  3  -  o)N\{o)do  jm 
=  JVr(-  +  /?  -  a). 

We  conclude  that 

Oi,r+i  :=  (o,.r  -  at+j,r  -r  ...  +  at+m-i,r)/m,  for  r  >  0.  (12.2) 

Here  is  the  full  algorithm. 

Algorithm  12.  Given  the  B-spline  coefficients  a  =  (a,)  of  s  €  Sj^zs,  its  B-spline 
coefficients  a  —  (a,)  with  respect  to  the  refined  knot  sequence  ZZ/m  can  be  computed  as 
follows: 

a"  j  ~  ® . m  —  1; 

m  —  1 

Q i,r  —  Qj- j,r-  1  / 1  r  —  2,...,k, 

J=0 

u ,  : “  (Lik¬ 
in  practice,  one  would  use  the  algorithm  repeatedly  with  m  -  2  rather  than  once  with 
a  larger  m.  For,  the  computational  cost  is 

nm(fc  -  l)((m  -  1 ) .4  -  D). 


with  n  the  number  of  coefficients  to  start  with,  and  A  and  D  the  cost  of  one  addition, 
respectively  division.  If,  e.g..  the  targeted  refinement  is  to  have  2Mn  coefficients,  then  the 


cost  ratio  of  the  choice  m  =  2M  versus  the  use  of  (j,  applications  of  the  algorithm,  each  time 
with  m  -  2,  is 


2"((2'i  -  1  )A  +  D) 

(2  +  2-  ~  ...  *  2»)[A  -  D) 


2 "''A  +  D/2. 


In  addition,  even  though  the  repeated  application,  with  m  —  2,  takes  roughly  twice  as 
many  divisions,  these  are  just  divisions  by  2. 


13.  Shape  preservation 

In  this  section,  we  use  knot  insertion  to  prove  the  shape  preserving  property  of  B- 
splines.  Roughly  speaking,  this  property  says  that  a  spline  has  the  same  shape  as  its 
control  polygon. 

We  begin  with  the 

convex  hull  property  If  tj  <  t  <  tj+u  then  $(t)  is  a  convex  combination  of  the  k 
B-spline  coefficients  ay_*+|, . . .  ,ar 

which  follows  from  Algorithm  9  or  directly  from  the  facts  that  B-splines  are  nonnegative 
(Sec.  1)  and  add  up  to  1  at  every  point  (see  (4.6)). 

For  a  statement  of  the  full  shape  preserving  property,  we  recall  that 

S' (a) 

is  the  standard  notation  for  the  number  of  (strong)  sign  changes  in  a  sequence  o.  Thus 
S"(l,  -  1, 1,-1)  =  3,  S"(l,0,l.-1)  =  1,  S~(0.0.0.0)  =  0. 

Theorem  13.  Variation  diminution  S  (s)  <  S-(fl)  :  i.e.,  with  xj  <  •  •  -  <  xT 
arbitrary, 

S~(s(x,),...,s(xr))  <  S~(a). 

Proof.  Recall  from  Sec.  11  that  s(xj), . . .  ,s(xr)  is  a  subsequence  of  the  sequence,  a 
of  coefficients  for  s  with  respect  to  the  refined  knot  sequence  t  which  contains  each  x,  at 
least  k  -  1  times.  Hence  it  is  sufficient  to  prove  that  S_(a)  <  S  (o).  But  this  follows  once 
we  know  that  S_(a)  <  S"(a),  with  a  obtained  by  (11.3).  i.e.,  by  insertion  of  just  one  knot. 
For  this  simple  case,  though,  the  conclusion  is  immediate  if  we  think  of  the  construction  of 
a  from  a  as  occurring  in  two  steps:  In  the  first  step,  we  insert  a,  between  a,  j  and  a,,  and 
this  does  not  increase  the  number  of  sign  changes  since  each  a,  is  a  convex  combination 
of  its  neighbors  a,_j  and  a ,  in  that  new  sequence.  In  the  second  step,  we  pull  out  a  as  a 
subsequence,  and  this  may  only  lower  the  number  of  sign  changes. 
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Corollary.  Shape  preservation  A  spline  crosses  any  straight  line  no  more  often 
than  does  its  control  polygon.  In  particular,  if  the  control  polygon  is  monotone  (convex), 
then  so  is  the  spline. 

Proof  Let  s  be  the  spline  and  p  the  straight  line.  Then  S~(s  -  p)  is  the  number 
of  times  the  spline  crosses  the  straight  line.  Since  s  -  p  is  a  spline,  this  is  bounded  by 
S~(a  -  6),  with  a,  6  the  B-spline  coefficients  of  s.  resp.  p  with  tespect  to  t,  and  this  equals 
the  number  of  times  the  control  polygon  Ca,t  crosses  the  control  polygon  for  p.  But,  as 
we  observed  in  Sec.  11.  the  control  polygon  for  the  straight  line  p  is  p  itself.  This  proves 
the  general  statement. 

For  the  particulars,  recall  that  a  (continuous)  function  is  monotone  if  and  only  if  it 
crosses  any  any  constant  function  at  most  once,  and  that  a  function  is  convex  if  it  crosses 
any  straight  line  at  most  twice  (dipping  first  below  and  then  rising  above  the  line  in  case 
it  crosses  it  twice). 


Figure  13.1  A  cubic  spline,  its  control  polygon,  and  various  straight  lines 
intersecting  them.  The  control  polygon  exaggerates  the  shape 
of  the  spline.  The  spline  crossings  are  bracketed  by  the  control 
polygon  crossings. 
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